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1.  INTRODUCTION 

The  Army  inventory  contains  several  models  which  compute  transmission  (T) 
through  an  obscuring  medium  composed,  for  example,  of  smoke  or  dust; ^  that  is. 


T  » 


e 


-T 


(1) 


where  t  is  the  optical  depth  along  the  path  of  propagation. 

On  one  hand  attempts  are  then  made  to  directly  relate  transmission  to  electro- 
optical  system  performance  and  smoke  effectiveness  by  considering  only  the 
directly  transmitted  signal: 


S{?)  =  Sir^lT  ,  (2) 

where  S(r)  is  the  optical  signal  received  by  an  observer  at  (r)  from  a  target 
at  fg.  The  transmission  (T)  includes  effects  of  both  scattering  out  of  the 
path  plus  absorption  along  the  path,  the  composite  process  being  referred  to 
as  extinction. 

On  the  other  hand  system  performance  modelers  know  that  electro-optical  sys¬ 
tems  (including  the  eye-brain)  respond  not  only  to  directly  transmitted  radia¬ 
tion  but  also  to  contrast,  the  definition  of  which  may  vary  among  models  but 
generally  requires  an  addition  to  equation  (2)  to  account  for  path  radiance, 
or  "brightness."  That  is  (see  also  figure  1), 

S(f)  .  S(f„)T  *  S|,(f)  .  ,3) 


where  the  contribution  due  to  path  radiance  (Sp)  may  be  due  either  to  scatter¬ 
ing  of  ambient  radiation  (for  example,  sun,  moon,  and  sky)  into  the  path  of 
propagation  or  (thermal)  emission  along  the  path,  or  both.  References  to 
scattering  out  of  and  into  are  emphasized  to  note  that  the  former  does  not 
directly  contribute  to  path  radiance  and  can  usually  be  treated  as  indistin¬ 
guishable  from  simple  Beer's  law  attenuation.  The  latter  however,  which  does 
contribute  to  path  radiance,  is  usually  a  complex  function  of  mapy  factors, 
including  angular  properties  of  both  the  scattering  medium  and  the  ambient 
radiation. 


iR.  A.  Sutherland,  D.  W.  Hoock,  and  R.  B.  Gomez,  1981,  An  Objective  Summary  of 
US  Amy  Electro-Optical  Modeling  and  Field  Testing  in  an  Obscuration  ^nvi ron- 
nieht~,  aSl-tR-009o,  05  Army  Atmospheric  Sciences  Laboratory,  White  Sands 
Missile  Range,  NM 
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Note  that,  unlike  transmission,  path  radiance  has  a  vector  nature  which  means 
physically  that  In  real  world  scenarios,  such  as  the  smoked  battlefield,  asym¬ 
metries  exist  between  target  and  observer,  giving  one  or  the  other  an  "optical 
advantage."  This  vector  nature  Is  the  essence  of  the  present  model  and  should 
not  be  overlooked  In  the  deceptively  simple  form  of  equation  (3)  or  by  the 
necessarily  complex  formulation  to  follow. 

The  existence  of  path  radiance  In  real  world  scenarios  Is  often  of  overriding 
significance  In  affecting  perception  and  Is  commonly  observed  In  nature.  One 
example  Is  the  apparent  disappearance  of  stars  In  daytime.  Another  Is  experi¬ 
enced  by  Individuals  driving  a  vehicle  through  fog  with  the  headlights  on  high 
beam.  In  both  cases  perception  Is  diminished  due  to  Interference  caused  by 
scattering  that  Is  manifested  by  path  radiance.  In  the  Infrared  the  effect  of 
path  radiance  Is  to  (partially)  offset  the  effects  of  absorption.  Another 
example  Is  radiance  data  sensed  via  orbiting  satellites.  Often  such  data  are 
highly  accurate  (~1®C)2  ^  when  Inverted  to  obtain  surface  temperature.  This 
accuracy  occurs  despite  the  fact  that  the  path  transmission  In  these  cases, 
even  In  the  so-called  atmospheric  "windows,"  Is  only  on  the  order  of  60  per¬ 
cent,  which  If  taken  alone  would  Imply  a  corresponding  temperature  error  on 
the  order  of  50  to  100®C!  The  explanation  here  lies  In  the  basic  physics  of 
Infrared  propagation  In  which,  for  practical  scenarios,  absorption  Is  always 
accompanied  by  Kirchhoff  (l.e.,  thermal)  emission  as  elucidated  for  the  case 
of  the  atmosphere  In  early  works.**  ^ 

The  degree  to  which  scattering  and/or  emission  can  be  Important  Is  Indicated 
by  the  optical  properties  of  the  medium;  the  best  Indicators  are  the  mass 
extinction  coefficient  (a)  which  Influences  total  extinction,  the  single 
scatterln^g  albedo  (uq)  which  Indicates  the  fractional  amount  of  scattering, 
and  1  -  wq  which  Indicates  the  fractional  amount  of  absorption. 


2r.  a.  Sutherland  et  al,  1979,  "A  Real  Time  Satellite  Data  Acquisition, 
Analysis  and  Display  System-A  Practical  Application  of  the  GOES  Network,"  J 
Appl  Meteorol ,  3:355-360 

^E.  Chen  et  al,  1979,  "Satellite-Sensed  Winter  Nocturnal  Temperature  Patterns 
of  the  Everglades  Agricultural  Area,"  J  Appl  Meteorol ,  8:992-1002 

**0.  D.  Kern,  1965,  "Evaluation  of  Infrared  Emission  of  Clouds  and  Ground  as 
Measured  by  Weather  Satellites,"  Environmental  Research  Papers,  No.  155, 
AFCRL-65-840,  Air  Force  Cambridge  Research  Laboratories,  Hanscom  Air  Force 
Base,  MA 

®S.  M.  Greenfield  and  W.  W.  Kellogg,  1960,  "Calculations  of  Atmospheric 
Infrared  Radiation  as  seen  from  a  Meteorological  Satellite,"  J  Meteorol, 
6:283-290  - 


Inventory  smokes  have  (dq+I  In  the  visible,^  indicating  a  predominance  of  scat¬ 
tering,  and  wq-K)  in  the  infrared,  indicating  a  predominance  of  absorption,  and 
consequently  emission.  Thus  path  radiance  is  important  and  perhaps  even  of 
overriding  significance  fo*"  inventory  smokes  from  the  visible  through  the 
infrared. 

The  need  for  further  model  development  in  this  area  was  established  in  an 
earlier  study^  which  made  a  detailed  examination  of  the  Amy  inventory  of 
existing  smoke  and  dust  obscuration  models.  A  major  finding  of  this  work  was 
that  although  most  models  reported  capabilities  for  treating  attenuation,  all 
were  deficient  in  wholly  treating  path  radiance  for  wavelengths  through  the 
infrared. 

As  a  step  toward  filling  this  technological  gap,  an  improved  smoke  obscuration 
model  reported  herein  was  developed.  Since  three  of  the  models  studied  (SOM 
II, ^  HECSOM,®  and  ACT-I,®  did  report  some  capabilities  in  the  visible,  the 
most  promising,  ACT- I,*  was  chosen  as  a  starting  point  (hence  the  acronym 
ACT-II  for  the  present  model). 

The  approach  is  to  provide  optical  information  critical  to  the  needs  of  pres¬ 
ently  existing  electro-optical  system  performance  and  smoke  effectiveness 
models.  An  infonTtal  survey  disclosed  that  the  requirements  were  reducible  to 
the  following  fundamental  quantities: 

1.  Ambient  irradiance  (light  level), 

2.  Target  and  background  radiance. 


®R.  C.  Shirkey  and  R.  A.  Sutherland,  1981,  "Aerosol  Phase  Function  Data  Base," 
chapter  16,  EOSAEL  80,  Volume  I,  Technical  Documentation,  editor  L.  D.  Duncan, 
ASL-TR-0072,  05  Army  Atmospheric  Sciences  Laboratory,  White  Sands  Missile 
Range,  NM  (AD  B055130L) 

^R.  A.  Sutherland,  D.  W.  Hoock,  and  R.  B.  Gomez,  1981,  An  Objective  Summary  of 
US  Army  El ectro-Opti cal  Modeling  and  Field  Testing  in~an  Obscurant  Environ¬ 
ment,'*  ASl-Tr-U096,  05  Army  Atmospheric  Sciences  Laboratory,  White  Sands 
Missile  Range,  NM 

^Smoke  Obscuration  Model  II  (SOM  11)  Computer  Code  Volume  II  -  Analyst  Manual, 
l9/9,  jTCG/ME  Smoke  and  Aerosol  working  Group  Document  61,  JYcG/ME-7&-9-Z 

®R.  K.  Dumbauld  and  H.  Bjorklund,  1977,  Mixing  Layer  Analysis  Routine  and 
Transport/Diffusion  Application  Routine  tor  tPAMS,  eCOM-zz-Z.  Atinospfieric 
Sciences  Laboratory,  Os  Amy  Electronics  Command,  White  Sands  Missile  Range, 
NM 

®R.  B.  Gomez,  R.  Pennsyle,  and  0.  Stadtlander,  1979,  "Battlefield  Obscuration 
Model,  ACT  I,"  Proceedings  of  Smoke  Symposium  III,  Harry  Diamond  Laboratories, 
Adel  phi,  MD  - - 


*The  acronym  ACT  derives  from  the  developing  agencies  Atmospheric  Sciences 
Laboratory,  Chemical  Systems  Laboratory,  and  TRASANA.  ~ 


3.  Line  of  sight  (LOS)  transmission,  and 

4.  LOS  path  radiance, 

where  target  and  background  radiances  are  computed  for  both  the  unperturbed 
(smoke  free)  and  smoked  environment,  and  LOS  data  are  provided  for  both  the 
observer-target  and  observer-background.  From  these  fundamental  quantities 
other  specialized  data  such  as  contrast  or  apparent  resolvable  temperature  can 
be  determined  easily  for  input  to  existing  system  performance  and  smoke  effec¬ 
tiveness  models  such  as  the  target  acquisition  model  of  the  Night  Vision  and 
Electro-Optics  Laboratory or  the  munition  expenditures  models  described  by 
Pennsyle^  and  Hoock.^^ 

In  another  respect  care  is  taken  so  that  the  (present)  model  inputs  are  com¬ 
patible  with  the  outputs  of  other  associated  models  such  as  LOWTRAN^^  and 
AGAUS^**  as  well  as  with  data  collected  during  field  tests  such  as  Smoke  Week 
11.15 

Although  the  primary  focus  is  on  optical  phenomena,  the  important  aspect  of 
obscurant  transport  and  diffusion  has  not  been  ignored.  The  approach  here  is 
to  generalize  procedures  so  that  the  model  will  accommodate  any  arbitrary 
ensemble  of  Gaussian  smoke  clouds,  providing  a  convenient  framework  for  possi¬ 
ble  future  union  with  equivalent  generalized  transport  and  diffusion  models. 
Most  present  models,  however,  do  not  provide  cloud  temperature  which  is  criti¬ 
cal  for  the  infrared.  Thus  a  Gaussian  diffusion  model  was  developed  based 


i°"Combat  Simulation  Target  Acquisition  Model  and  Data  Input"  (U), 
CONFIDENTIAL,  1980,  Draft  Technical  Report,  US  Army  Night  Vision  and  Electro- 
Optics  Laboratory,  Fort  Belvoir,  VA  (in  process) 

Hr.  0.  Pennsyle,  1979,  Methodol o^  Smoke/Obscurant  Muni ti on 
Expenditure  Requirements'!  aRCsl-TR-7902z,  Chemical  Systems  Laboratory, 
Aberdeen  Proving  Ground,  MD 

Ho.  W.  Hoock,  1981,  "SCREEN,"  chapter  5,  EOSAEL  80,  Volume  1,  Technical 
Documentation,  editor  L.  D.  Duncan,  ASL-TR-007Z,  uS  Army  Atmospheric  sciences 
laboratory.  White  Sands  Missile  Range,  NM  (AD  B055130L) 

Hj.  E.  Selby  et  al,  1978,  "Atmospheric  Transmittance/Radiance:  Computer  Code 
LOWTRAN  4,"  Environmental  Research  Papers,  No.  626,  AFGL-TR-78-0053,  Air  Force 
Geophysics  Laboratory,  Hanscom  Air  Force  Base,  MA 

^“♦R.  C.  Shirkey  et  al,  1980,  Single  Scattering  Code  AGAUSX:  Theory,  Applica¬ 
tions,  Comparisons,  and  Listing,  asl-TR-006z,  Us  Army  Atmospheric  Sciences 
Laboratory,  white  Sands  Missile  Range,  NM 

Hqpg  Final  Test  Report  on  Smoke  Week  II  at  Eg! in  AFB,  FL  (U),  CONFIDENTIAL, 
1978,  Volumes  1  and  II,  DPG-Fr-76-317,  Dugway  Proving  Ground,  UT 


upon  commonly  used  procedures'^  and  extended  to  include  buoyant  rise  and 
cloud  temperature  using  fundamental  principles.^® 

The  work  is  divided  into  three  parts:  the  present  work  covers  theory  and 
examples,  a  second^o  covers  program  documentation  and  a  users  guide,  and  a 
third^i  covers  validation  and  applications. 

Historically  the  problem  of  path  radiance  and  its  significance  to  visible 
perception  have  been  recognized  by  the  Army  modeling  community  for  several 
years.  As  early  as  1972  an  unpublished  document  described  a  smoke  obscuration 
model  (SOM)  which  reported  to  compute  visible  contrast  and  was  later  accepted 
as  the  Joint  Technical  Coordinating  Group  (JTCG)  working  model.  This  early 
model  was  expanded  by  at  least  two  groups,  one  leading  to  the  development  of 
the  model  SOM  and  another  to  ASLSOM  which  was  further  modified  to  become 
the  ACT  model®  which  is  the  direct  forerunner  of  the  present  model. 


2.  OUTLINE  AND  SCOPE 

The  fundamental  optical  quantities  to  be  determined  in  addition  to  transmis¬ 
sion  are  the  amounts  of  radiant  energy  received  by  an  observer  from  the  two 
directions  (approximately  coincident)  defined  by  the  relative  positions  of  a 


^®Smoke  Effectiveness  Manual,  1979,  JTCG/ME  Smoke  and  Aerosol  Working  Group 
Document  Number  FM  101-61-8 

^'^F.  V.  Hansen,  1979,  Engineering  Estimates  for  the  Calculation  of  Atmospheric 
Dispersion  Coefficients,  ASL  Internal  Report,  US  Army  Atmospheric  Sciences 
Laboratory,  White  Sands  Missile  Range,  NM 

^®F.  Pasquill,  1974,  Atmospheric  Diffusion,  second  edition,  Halsted  Press 
Div.,  John  Wiley  and  Sons,  Inc.,  New  York 

^®C.  H.  B.  Priestley,  1956,  "A  Working  Theory  of  the  Bent-Over  Plume  of  Hot 
Gas,”  Quart  J  Roy  Meteorol  Soc,  82:165-176 

2®R.  a.  Sutherland  and  D.  Clayton,  1981,  An  Improved  Smoke  Obscuration  ^del 
Act  II:  Part  2  Documentation  and  User  Guide,  Technical  Report,  US  Arnty  Atmo¬ 
spheric  Sciences  Laboratory,  White  Sands  Missile  Range,  NM  (in  process) 

21r.  a.  Sutherland,  1981,  "Comparisons  Between  the  Improved  Smoke  Obscuration 
Model  ACT  II  and  Recent  Smoke  Week  Data,"  Proceedings  of  Smoke  Symposium  V, 
Harry  Diamond  Laboratories,  Adel  phi,  MD 

^Smoke  Obscuration  Model  II  (SOM  II)  Computer  Code  Volume  II  -  Analyst  Manual , 
1979,  jTCG/ME  Smoke  and  Aerosol  Working  Group  Document  6l,  jYcg/mE-78-9-2 

®R.  B.  Gomez,  R.  Pennsyle,  and  D.  Stadtlander,  1979,  "Battlefield  Obscuration 
Model,  ACT  I,"  Proceedings  of  Smoke  Symposium  III,  Harry  Diamond  Laboratories, 
Adel  phi,  MD 
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target  and  background,  both  treated  as  Lambertian  surfaces.  The  radiance 
Incident  at  the  observer  from  each  direction  Is  generally  composed  of  two 
parts:  (1)  the  direct  radiance  emitted  and  reflected  by  the  target  (or  back¬ 
ground)  then  transmitted  (with  some  loss  due  to  extinction)  along  the  LOS  to 
the  observer  and  (2)  the  diffuse,  or  path,  radiance  emitted  and  scattered  by 
suspended  material  (smoke)  at  ai)  points  (such  as  P  In  figure  2)  along  the  LOS 
then  transmitted  (again  with  some  loss  due  to  extinction)  a  remaining  distance 
to  the  observer,  giving  rise  to  a  path  radiance.  One  aspect  of  the  problem 
which  causes  major  complexity  Is  that  the  entire  environmental  sphere  must  be 
considered  the  source  for  the  reflected  and  diffuse  radiation,  thus  requiring 
Integration  over  all  angles,  not  only  at  the  target  and  background  but  also  at 
all  points  along  the  LOS.  Except  for  rare  circumstances,  these  Integrations 
must  be  carried  out  by  some  approximate  numerical  technique.  The  approach 
taken  here  Is  to  divide  the  sky  hemisphere  Into  discrete  angular  sectors  and 
then  assume  that  the  radiances  from  the  various  sectors  are  either  known  from 
measurement  (as  In  recent  field  tests)  or  produced  by  some  appropriate  model 
(perhaps  LOWTRAN^^).  Terrain  radiance  due  to  reflected  sky  radiation  and 
thermal  emission  can  then  be  calculated  from  knowledge  of  surface  albedo, 
emissivlty,  and  temperature  to  complete  the  characterization  of  the  (smoke 
free)  radiation  environment,  which  is  then  assumed  constant  throughout.  The 
exact  sectoring  procedure  used  In  the  model  Is  outlined  In  section  3. 

Mathematically  the  problem  can  be  summarized  by  the  following  formal  expres¬ 
sion  describing  radiant  propagation  along  a  straight  path  over  a  distance  r.22 


R(?) 


R(ro)e"^^’"’ 


r 

I 


C0(r' ,  e,  (f)) 


(1  -  2;o)B(x,  Tp)] 


,-T(r, 


r') 


dr’ 
Hr^  O'" 


(4) 


where  R(r)  is  the  radiance  Incident  at  r;  Rir^)  is  the  radiance  of  the  target 
(or  background)  located  at  fg;  and  9,  <j)  are  the  polar  angles  defining  the  path 
of  propagation  (that  Is,  the  LOS).  Although  the  vector  notation  will  be 
dropped,  it  Is  assumed  here  and  In  the  following  that  the  observer  is  at  the 
origin  and  that  the  coordinates  are  rotated  so  that  r  and  f'  lie  along  the  LOS 
(figure  4).  Generally  the  term  J(r,  9,  (ji)  Jcalled  the  source  function) 
accounts  for  scattering  Into  the  LOS,  and  (1  -  uo)B(x,  T  )  accounts  for  emis¬ 
sion  from  Increments  along  the  LOS.  ^ 


E.  Selby  et  al,  1978,  "Atmospheric  Transmittance/Radiance:  Computer  Code 
LOWTRAN  4,"  Environmental  Research  Papers,  No.  626,  AFGL-TR-78-0053,  Air  Force 
Geophysics  Laboratory j  Hanscom  Air  Force  Base,  MA 

2^S.  Chandrasekhar,  1960,  Radiative  Transfer,  second  edition,  Dover  Press 
Publications,  Inc.,  New  York 
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The  Planck  or  blackbody  function  of  equation  (4)  Is  written  explicitly  as 


B(x,  Tp)  = 


2hc^\~^  AX 
[exp{hc/kTp)  -  1] 


(5) 


where  x  Is  the  wavelength,  ax  the  bandpass,  and  h,  k,  and  c  are,  respectively, 
the  Planck  constant,  the  Boltzmann  constant,  and  the  speed  of  light  In  vac¬ 
uum.  The  obscurant  temperature  (Tp)  Is  assumed  variable  over  the  path,  so 
that  B  contains  an  Implicit  dependence  on  r. 

The  optical  thickness  (t)  Is  defined  as 


T(r.  r') 


oC(r")dr'' 


(6) 


where  o  is  the  obscurant  mass  extinction  coefficient  and  C  is  the  obscurant 
concentration.  Both  obscurant  temperature  and  concentration  are  discussed  in 
section  5. 

The  source  function  J{r,  0,  <>)  Is  difficult  to  compute,  requiring  Integrations 
over  the  entire  environmental  sphere  accounting  for  the  angular  characteris¬ 
tics  of  both  the  ambient  radiation  and  the  scattering  medium.  Except  for 
trivial  cases,  no  exact  methods  exist  for  computing  this  term;  and  for  real¬ 
istic  scenarios,  some  approximate  technique  must  be  employed.  The  model  uses 
the  single  scattering  approximation  In  which  the  source  function  can  be 
written  as 


0(r,  9,  t) 


/  P(0J  L(0',  ’'s'dn' 

4ir  ^ 


(7) 


where  e.  Is  the  scattering  angle  (figure  2)  and  P  the  phase  function.  The 
term  L(0',  (|)’)  consists  of  two  parts:  the  source  radiance  Rj(0',  (p')  from 
the  directions  of  the  sky  and  terrain  sectors  and  the  thermal  emission  along 
these  same  directions.  Mathematically, 


L{0',  i>')  =  R5(0',  ((.')  +  (1  -  uq) 


B(X,  Tp)e^^*'"’  ’"'^^dr' 


(8) 


In  the  above  expressions,  r  Is  distance  to  aqy  point  along  the  LOS;  r'  Is 
distance  from  that  point  along  the  direction  defined  by  0',  4';  r^  Is  distance 
to  the  sky  and  terrain  sources;  and  dn'  is  the  differential  solid  angle. 


For  Inventory  smokes  (and  neglecting  polarization),  the  angular  dependence  of 
the  phase  function  Is  dependent  only  upon  the  scattering  angle,  which  from 
simple  geometry  Is  given  by  (see  figure  2): 


cos  =  [cos  0  cos  0'  +  sin  0  sin  0'  cos((|)  -  (!>')] 


(9) 


Some  caution  Is  required  in  using  equation  (9)  to  assure  the  proper  algebraic 
sign.  For  use  In  the  phase  function  equation  (9)  is  correct  as  it  stands,  but 
for  Lambertian  surfaces  (that  is,  target,  background,  etc.)  the  sign  must  be 
reversed  because  the  convention  used  in  the  model  requires  the  surface  normal 
pointing  positive  inward  (for  example,  away  from  the  observer)  which  in  turn 
requires  the  reversal  in  sign. 

The  phase  function  is  required  as  Input  but  can  readily  be  obtained  from  the 
associated  model  AGAUS,^**  one  version  of  which  is  distributed  with  the 
Electro-Optical  Systems  Atmospheric  Effects  Library  (EOSAEL  80^^).  The  phase 
function  is  assumed  to  be  normalized  such  that 


J  P(0^)dn  =  u)o  , 

4ir 


but  in  the  model  it  is  renormalized  via  equation  (10)  to  a  single  scattering 
albedo  specified  as  input.  However,  to  be  strictly  compatible  with  theory  the 
input  single  scattering  albedo  should  be  that  computed  from  Mie  scattering. 

The  major  objective  of  the  model  is  to  evaluate  the  two  components  of  equation 
(4),  once  for  the  observer-target  and  once  for  the  observer-background  by 
using  the  procedures  described  by  equations  (4)  through  (9).  For  the  special 
case  of  computing  R(rQ),  the  target  or  background  radiance,  the  same  procedure 
for  the  second  term  of  equation  (4)  is  used  except  that  the  factor  ^ 
equation  (7)  is  replaced  by  (a  cos  Q^/u)  which  assumes  a  Lambertian  surface  of 
albedo*  (a)  with  surface  normal  along  the  LOS.  Also  for  these  cases  the  com¬ 
putations  are  restricted  to  0^  ^  90®  to  avoid  contributions  due  to  reflection 
from  the  rear  surface.  The  (smoke  free)  surface  irradiance  (^sfc)  is  also 
computed  in  the  same  manner  with  0  =  (^  =  0  (vertical)  and  the  factor  lA 


^‘♦R.  C.  Shirkey  et  al,  1980,  Single  Scattering  Code  AGAUSX:  Theory,  Applica- 
t1 ons ,  Compari sons ,  and  Li sti ng ,  ASL-tR-ODfeZ ,  US  Army  Atmospheric  Sciences 
Laboratory,  white  Sands  Missile  Range,  NM 

23r.  a.  Sutherland,  1981,  "Smoke  Obscuration  Model,"  chapter  3,  EOSAEL  80, 
Yol ume  1 ,  Techni cal  Documentati on,  editor  L.  D.  Duncan,  ASL-TR-00/2,  us  Army 
Atmospheric  Sciences  Laboratory,  White  Sands  Missile  Range,  NM  (AD  B055130L) 

♦For  opaque  surfaces,  reflectivity  (r),  albedo  (a),  and  emissivity  (c)  are 
related  as  (a  =  r  and  c  =  1  -  r) . 
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removed.  In  all  cases,  an  emission  term  of  the  form  eB(x,  T)  is  added  where  T 
is  chosen  appropriately  as  the  surface,  target,  or  background  temperature,  and 
the  emissivity  (e)  is  computed  from  the  reflectivity  or  albedo  as  e  =  (1  - 
a).  The  (smoke  free)  surface  irradiance  is  used  later  (see  equation  (18))  to 
compute  radiances  for  terrain  sectors  which  are  then  treated  in  the  same 

manner  as  sky  sectors. 

For  the  visible  scenarios,  the  effect  of  emission  will  be  negligible  because 
of  the  smallness  of  the  blackbody  function  in  these  spectral  regions  at 

nominal  temperatures.  For  infrared  scenarios,  this  term  often  dominates, 
being  more  pronounced  at  higher  temperature,  which  means  that  errors  due  to 
neglect  of  multiple  scattering  will  be  minimal  in  the  infrared.  However, 
errors  may  occur  in  the  infrared  due  to  uncertainties  in  the  cloud  tempera¬ 
ture. 

The  process  to  be  modeled  here  can  be  summarized  in  geometrical  terms  with  the 
aid  of  figure  2.  Simply  stated,  the  problem  is  to  compute  contributions  to 
path  radiance  at  each  point  P  along  the  LOS,  and  then  to  sum  over  all  such 

points.  At  each  increment,  effects  of  extinction  must  be  included  over  the 

remaining  path  TO  to  the  observer.  At  each  point  P  the  contribution  is  com¬ 
posed  of  two  parts--one  due  to  scattering  into  the  increment  from  all  angles 
and  the  other  due  to  emission  by  the  increment.  The  single  scattering 
approximation  assumes  that  the  radiance  along  any  path  "SP  is  scattered  into 
the  LOS  only  once  and  that  this  scattering  occurs  at  P.  Thus  the  radiance 
scattered  into  the  LOS  at  point  P  consists  of  the  source  radiance,  R^,  reduced 
by  extinction  over  the  path  plus  the  summation  of  the  emission  from  each 
element  P'  along  the  roth  the  emissive  contribution  of  each  element  is 
reduced  by  extinction  over  the  path  PP'.  The  total  scattering  contribution  of 
each  increment  at  P  is  found  by  summing  over  all  angles,  accounting  for 
angular  scattering  properties  of  the  medium  via  the  phase  function.  Total 
path  radiance  is  found  by  summing  over  all  increments  along  OT. 

In  the  model  the  increment  spacings  are  chosen  by  a  criterion  based  upon  the 
incremental  optical  depth  at.  This  method  speeds  computations  by  avoiding 
insignificant  contributions  for  increments  containing  no  obscurant  which  would 
occur  for  a  criterion  based  on  spatial  separation  (Ar).  The  minimum  spacing 
in  the  model  however  is  normally  defaulted  to  1  m. 

The  model  treats  extinction  due  to  the  ambient  atmosphere  by  appropriately 
modifying  transmission  (i.e.,  T^qj  =  T^^Qj^g  ^atmos^  along  the 

LOS.  This  option  is  employed  by  way  of  a  user  supplied  volume  extinction 
coefficient,  (a),  so  that  T^^g^  *  where  L  is  distance  of  propagation. 

Parallel  point  sources  of  radiation,  including  the  sun  or  moon,  are  also 
treated  by  the  model . 

In  all  of  the  above  computations,  the  model  computes  optical  thickness  (t)  by 
assuming  the  medium  composed  of  any  ensemble  of  obscuring  smoke  clouds  defined 
by  centroid  locations,  Gaussian  standard  deviations,  and  temperature.  Methods 
for  integrating  equation  (6)  and  for  producing  the  ensemble  are  given  in  later 
sections  of  the  report. 
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3.  AMBIENT  IRRAOIANCE  (SECTORING  SCHEME) 


This  section  describes  the  sectoring  scheme  used  to  simulate  incoming  radia¬ 
tion  from  sky  and  terrain  which  will  then  be  used  to  approximate  terms  for  the 
source  function  of  equation  (8).  Throughout  this  section  repeated  use  will  be 
made  of  approximations,  assuming  that  scenario  relative  distances  are  small  In 
comparison  to  spatial  variations  In  ambient  conditions.  This  process  consid¬ 
erably  simplifies  the  geometry  by  allowing  all  scenario  elements  to  be  treated 
as  exposed  to  the  same  ambient  radlatlonal  environment.  These  are  approxima¬ 
tions  often  used  In  problems  of  this  type  and  Introduce  only  minimal  errors. 

The  major  divisions  of  the  entire  4ir  steradians  comprising  the  environmental 
sphere  are  sketched  In  figure  3.  The  upper  sector  Is  assumed  to  be  comprised 
of  sky  (Including  sun,  moon,  and  clouds)  and  the  lower  to  be  overall  flat  ter¬ 
rain.  Both  sky  and  terrain  will  be  treated  as  sources  of  ambient  radiation, 
the  latter  through  reflection  of  sky  radiation  and  thermal  emission. 

To  facilitate  computations,  the  two  major  regimes  are  further  subdivided  Into 
angular  sectors  subtending  equal  solid  angles.  These  discrete  sectors  are 
then  treated  as  point  sources  of  parallel  radiation  emanating  from  the  direc¬ 
tion  of  the  sector  midpoint.  Additional  sources  of  radiation  such  as  the  sun 
or  moon  are  superimposed  at  their  appropriate  angular  positions.  The  model 
will  accommodate  variable  radiance  from  each  of  the  discrete  sky  sectors,  but 
to  maintain  consistency  with  the  assumptions  mentioned  earlier,  one  must 
assume  that  the  terrain  Is  homogeneous  In  albedo,  emissivlty,  and  temperature. 

The  procedure  for  sectoring  the  two  regimes  Into  equal  angular  sectors  follows 
directly  from  the  definition  of  solid  angle;  dn  =  sin  e  de  d(}>,  where  e  and 
4i  are  the  usual  zenith  and  azimuth  angles.  ’ 

The  azimuthal  sectoring  is  particularly  simple  since  Integration  over  contig¬ 
uous  divisions  (((i^,  4'^  +  ^)  yields,  simply,  dn^  =  A())  sin  9  de,  where  = 

'*’1+1  "  '*’1  azimuthal  separation,  which  for  m  sectors  Is  simply  A4)  = 
2Tr/m.  The  representative  midpoints  are  then 


=  (1  -  l)iT/m  . 


(11) 


For  the  zenith  sectors,  the  Integration  between  contiguous  divisions  yields 

An  =  A(ti(C0S  9j  -  cos  Oj  +  j)  •  (12) 

For  n  sectors,  all  of  which  are  equal  and  contained  In  a  total  solid  angle  2n, 
we  have 


A((>(cos  0j  -  cos  Qj+j)  *  Zir/nm 


(13) 
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which  after  substituting  for  and  rearranging  becomes 


cos  *  cos  Oj  -  1/n  , 

from  which  all  divisions  can  be  calculated  by  knowing  that  Oj  =  0. 
lent  but  sometimes  more  convenient  expression  Is 

cos  0^  =  1  -  (j  -  l)/n  . 

J 


(14) 
An  equi va¬ 
ns) 


Further  reasoning  yields  the  following  equation  for  sector  midpoints; 

cos  "Sj  =  1  -  {2j  -  l)/2n  . 


(16) 


The  corresponding  distances  to  the  terrain  sector  midpoints  are 


r^’*^  =  h/cos  0  . 

^  J 


(17) 


where  h  Is  the  vertical  distance  from  the  surface  for  the  particular  scenario 
element  under  consideration.  The  radiance  from  the  sector,  assuming  a 
Lambertian  surface  Is 


R  =  [(a/Tr)E^^^  +  (1  -  a)B(x,  Tg^^)] 


(18) 


where  a,  and  T^^^.  are,  respectively,  the  surface  albedo,  Irradlance  (see 

section  2),  and  temperature.  From  this  point  on,  the  only  difference  In 
treating  si^  or  terrain  sectors  Is  that  the  finite  distance  to  the  terrain 
sectors  must  be  considered  via  equation  (17),  whereas  the  sky  sectors  can  be 
assumed  at  infinity  (actually  10,000  m  In  the  model). 

Ordinarily  one  does  not  have  sufficient  data,  or  the  inclination,  to  provide 
the  radiance  values  for  all  of  the  sectors  used  In  the  model;  therefore,  the 
model  was  programmed  to  proportion  the  sectors  uniformly  by  Interpolating  the 
Input  radiance  values  from  arbitrary  angles.  This  Interpolation  makes  the 
model  Input  directly  compatible  with  sky  radiance  data  from  the  smoke  tests. 


r 


Also,  to  avoid  Inconsistencies  between  the  computed  surface  Irradlance 

and  the  reported  measurements,^^  the  sun  and  sky  Input  data  are  treated  only 
as  relative  and  are  normalized  so  as  to  reproduce  the  measured  value  when 
Integrated  over  the  sky  hemisphere.  Thus  the  model  as  now  coded  requires  only 
relative  data  from  sun  and  sky  but  an  absol ute  determination  of  surface  Irra¬ 
dlance.  In  effect  this  method  reduces  the  complexity  of  the  required  Input. 


4.  OPTICAL  THICKNESS  CALCULATIONS 

This  section  describes  the  general  method  used  to  compute  smoke  concentration 
C(r)  and  optical  thickness.  Throughout  we  will  assume  a  constant  extinction 
coefficient  so  that  the  optical  depth  is  simply  the  product  (aCL)  where  CL  is 
the  line  Integrated  concentration,  commonly  called  CL  product. 

The  methodology  is  based  upon  the  general  assumption  that  a  smoke  plume  or 
cluster  can  be  represented  by  a  series  of  spatially  and  temporally  discrete 
overlapping  clouds  each  with  concentration  given  by  a  tri variant  Gaussian 
function.  This  Is  a  common  assumption  used  In  many  models  although  the  manner 
of  spacing  and  sizing  such  clouds  may  vary  from  model  to  model.  For  this 
latter  reason  the  methodology  is  kept  general  so  as  to  be  easily  adaptable  to 
various  cloud  transport  models. 

For  some  i^*’  cloud  centered  at  7. ,  T. ,  and  7.  the  concentration  (due  to  this 
source  only)  Is  given  by 


C^(x,  y,  z)  =  (2ir)"  ^^2  — I — ^  e"  ^2 

®x  ®y  ‘'z 

(19) 


where  Is  the  total  mass  of  the  cloud  and  accounts  for  (1)  munition  fill 
mass  expended  during  the  burn  producing  the  cloud,  (2)  munition  efficiency, 
and  (3)  smoke  yield  factor.  The  total  concentration  Is  found  by  summing  the 
concentrations  of  all  such  clouds. 


21r.  a.  Sutherland,  1981,  "Comparisons  Between  the  Improved  Smoke  Obscuration 
Model  ACT  II  and  Recent  Smoke  Week  Data,"  Proceedings  of  Smoke  Symposium  V, 
Harry  Diamond  Laboratories,  Adel  phi,  MD 
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It  is  convenient  to  rewrite  equation  (19)  in  spherical  coordinates  to  give  an 
expression  for  concentration  along  an  LOS  defined  as  before  by  polar  angles  e 
and  t  at  some  arbitrary  point  a  distance  r  from  the  origin.  It  is  straight¬ 
forward  to  show  that  the  equivalent  to  equation  (19)  is 


C(r,  0,  = 


a  a  a 
X  y  2 


(=e)'  •  (¥1  • 


where  the  indices  have  been  dropped  to  avoid  cumbersome  notation.  The  line 
parameters  a,  6,  y  and  offsets  x^,  y^,  are: 


o  =  sin  9  sin  41 


Xq  =  ^i  -  xi 


e  =  sin  0  cos  41 


r  =  cos  0 


yo  =  7,  -  y, 
^ 


where  xi,  yi,  and  21  are  coordinates  of  any  point  on  the  LOS,  taken  in  the 
model  to  be  the  common  point  such  as  P  in  figure  4. 

With  considerably  more  algebraic  manipulation  which  involves  expanding  the 
expression  in  the  exponential,  rearranging  and  then  rewriting  the  resultant 
expression  as  a  perfect  square,  the  following  expression  results: 


C^(r,  9,  4*)  =  Q|  e' 


1/.  r  -  R. 


which  is  itself  a  Gaussian  with  mean  standard  deviation  and  strength 
given  by  the  following  expressions,  again  with  indices  suppressed: 


(23)  cont 


The  final  desired  result  for  line  integrated  concentration  beginning  at  point 
P(xi,  yi,  zi)  along  the  line  described  by  (a,  e,  r)  for  a  distance  D  becomes 
(see  figure  4): 


CL^.(D,  0, 


(24) 


where  the  error  function  is  defined  as 

X 

erf(x)  =  ^  exp(-t2)dt  125) 

0 

and  is  computed  in  the  model  according  to  the  approximate  technique  as 
described  by  Abramowitz  and  Stegun.^** 

The  formulation  here  applies  in  a  wind  vector  aligned  coordinate  system 
requiring  that  scenario  Cartesian  coordinates  and  angles  be  transformed  to 
this  system  before  the  calculations. 

The  model  assumes  both  concentration  and  temperature  to  be  Gaussian  so  that  a 
relationship  analogous  to  equation  (24)  is  used  to  obtain  temperature  of 
various  line  segments  for  computation  of  thermal  emission.  Also  symmetric 
"image"  clouds  accounting  for  surface  particulate  reflection  are  included  in 
the  usual  manner. ^3 


5.  CLOUD  CONCENTRATION  AND  TEMPERATURE 

The  preceding  sections  assumed  a  transport  and  diffusion  model  generating  some 
pattern  of  overlapping  Gaussian  clouds.  Several  methodologies  which  can  be 


^‘♦M.  Abramowitz  and  I.  Stegun,  1970,  Handbook  of  Mathematical  Functions,  Dover 
Press  Publications,  Inc.  New  York 

23R.  a.  Sutherland,  1981,  "Smoke  Obscuration  Model,"  chapter  3,  EOSAEL  80, 
Volume  1,  Technical  Documentation,  editor  L.  D.  Duncan,  ASL-TR-007Z,  uS  Army 
Atmospheric  Sciences  Laboratory,  White  Sands  Missile  Range,  NM  (AD  BOSSISOL) 
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adapted  to  this  general  concept  are  available.  We  borrow  bits  and  pieces  from 
these  methodologies  to  produce  a  submodel  to  be  used  for  the  validation 
studies  reported  later.  Production  of  this  model  consists  of  generating  the 

parameters  Q^,  (IT^,  7^,  (o^,  o^,  o^),  and  cloud  temperature  which  will 

now  be  covered  In  order. 


5.1  Smoke  Source  Function  (Q) 

The  factor  Q  represents  the  total  mass  of  a  smoke  cloud  and  Is  composed  of  the 
product  of  factors  M,  x,  and  Y  where  M  is  the  mass  of  munition  fill  expended 
during  the  burn  producing  the  cloud*  X  is  the  chemical  efficiency  with  which 
the  mass  Is  converted  to  actual  smoke  nuclei,  and  Y  Is  the  smoke  yield  factor 
which  accounts  for  Increased  mass  due  to  hygroscopic  Interactions  with  the 
ambient  air  mass. 

For  Instantaneous  bursts  such  as  bulk  fill  white  phosphorus  munitions,  a 
single  cloud  of  mass  Q  =  M  X  Y  Is  used.  For  munitions  of  extended  burning 
time  (>  Is),  the  plume  Is  generated  as  a  series  of  discrete  puffs  produced 
during  short  time  Increments  (nominally  1  s).  Variable  burn  rate  Is  Included 
by  employing  either  a  quadratic  or  exponential  function  with  coefficients  as 
determined  empirically  from  field  tests.  The  EOSAEL  80  Technical  Documenta- 
t1on23  contains  a  review  of  these  burn  coefficients  and  other  munition  charac¬ 
teristics. 


5.2  Cloud  Centroids  (7,  T,  and  D 

With  the  coordinate  system  rotated  to  align  the  positive  x  axis  along  the  wind 
vector  and  assuming  the  cloud  to  be  transported  by  the  mean  wind  (O'),  the 
cloud  centroids  are  modeled  as 


^  “  ^m  (26) 

I  =  +  H(t) 

where  Y^u,  and  are  munition  coordinates.  The  method  of  computing  the 

cloud  rise  function  H(t)  which  also  Involves  the  cloud  temperature  Is  dis¬ 
cussed  later.  The  mean  windspeed  Is  computed  by  averaging  vertically  over  the 


23r.  a.  Sutherland,  1980,  "Smoke  Obscuration  Model,"  chapter  3,  EOSAEL  80, 
Volume  1,  Technical  Documentation,  editor  L.  D.  Duncan,  ASL-TR-0072,  US  Army 
Atmospheric  Sciences  Laboratory,  White  Sands  Missile  Range,  NM  (AD  B055130L) 


21 


significant  cloud  extent  (3a)  using  the  usual  vrindspeed  power  law:  U(z)  = 
Uj.(z/Zj.)^  where  Is  the  windspeed  at  an  (Input)  reference  height  (Z^)  and  P 
Is  the  vertical  profile  exponent. 


5.3  Dispersion  Functions  (a  ,  o  ,  a  ) 

A  y  4 

The  dispersion  functions  o  ,  o  ,  a  are  all  expressed  as  power  functions  of 

X  y  z 

the  X  centroid  with  Initial  offset,  a(0);  that  Is, 

°X.  y.  z  ° 

The  source  sigmas  a(0),  essentially  representing  the  dimensions  of  the  cloud 
at  t  =  0  are  modeled  by  the  following  power  functions  which  were  derived  from 
the  data  of  AMSAA  TR-201.25 


■ 


(7  in  meters) 


(27) 


0.3 

JO)  =  5.0  Q 
^ »  y 


0  .3 

0^(0)  =  1.7  Q 


(Q  in  kilograms) 


(28) 


The  diffusion  parameters  A  and  B  of  equation  (27)  are  modeled  as  functions  of 
the  surface  average  roughness  element  (Z^)  and  the  stability  category  as 

listed  in  table  1.  For  surface  roughness  Zq  >  0,  the  values  are  those  cited 

by  Hansen and  for  these  cases  the  parameter  C  of  equation  (27)  is  set  to 

zero.  For  a  roughness  entered  as  Z^  7  0  (default),  the  method  of  the  Smoke 

Effectiveness  Manual^®  is  used,  in  which  case  the  initial  sigmas  are  absorbed 

in  the  parameter  C,  and  the  term  a(0)  is  set  to  zero. 


Analysis  of  the  Smoke  Cloud  Data  from  the  August  1975  Jefferson  Proving 
Ground  Smoke  Test,  WIT,  AMSAA  Technical  Report  tR-zOl,  Aberdeen  Proving 
Ground,  MB'  CM)"  TO45874) 

^^F.  V.  Hansen,  1979,  Engineering  Estimates  for  the  Calculation  of  Atmospheric 
Dispersion  Coefficients,  ASL  Internal  Report,  US  Army  Atmospheric  Sciences 
Laboratory,  White  Sands  Missile  Range,  NM 

Smoke  Effectiveness  ^nual ,  1979,  JTCG/ME  Smoke  and  Aerosol  Working  Group 
Document  Humber  KM  lDl-61-8 
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r 


Following  the  methodology  cited  by  Hansen, a  ,  o  ,  and  a  are  reduced  by 

A  y  £ 

factors  0.74,  0.67,  and  0.67,  respectively,  for  Instantaneous  sources. 


5.4  Cloud  Temperature  and  Buoyant  Rise 


Current  methods  for  modeling  buoyant  rise  are  generally  limited  to  empirical 
methods  based  upon  observations  of  factory  smoke  stack  effluents or  curve 
fits  to  data  from  field  tests. These  procedures,  although  of  approximate 
validity  for  special  circumstances,  have  severe  shortcomings  for  the  general 
case  where  it  becomes  necessary  to  simultaneously  model  cloud  temperature 
consistently.  This  consistency  Is  particularly  important  In  the  infrared 
where  cloud  temperature  acquires  an  added  significance  of  its  own  1n  addition 
to  the  Indirect  effect  on  buoyancy.  The  method  developed  for  the  model 
applies  basic  principles  and  certain  simplifying  assumptions  borrowing  heavily 
from  earlier  works in  a  self-consistent  numerical  scheme  as  outlined 
below. 


The  buoyant  motion  Is  modeled  by  treating  each  cloud  of  the  ensemble  as  though 
Independent  of  other  clouds,  an  assumption  consistent  with  the  transport  and 
diffusion  methodology  discussed  earlier.  Initial  cloud  temperature  Is  modeled 
by  equating  the  Internal  thermal  energy  of  each  instantaneous  cloud  to  the 
energy  expended  during  the  exothermal  reaction  producing  the  cloud.  Assuming, 
as  before,  similar  distributions  In  both  temperature  and  concentration,  the 
following  expression  results  for  initial  cloud  temperature: 


T 


c 


(29) 


where  E  Is  the  obscurant  heat  of  reaction  (calorie/gram),  the  mean  concen¬ 
tration,  pCp  the  volumetric  specific  heat  of  the  ambient  air  (290  cal  m"^ 
"C"M,  and  T^  the  ambient  air  temperature  at  the  cloud  centroid.  The  use  of 
equation  (29)  assumes  complete  thermal  mixing  between  cloud  and  entrained  air. 


A.  Briggs,  1965,  "A  Plume  Rise  Model  Compared  with  Observations,"  J  Air 
Poll  Control  Assoc,  15:433 

^^Analysls  of  the  Smoke  Cloud  Data  from  the  1975  Je^fferson  Proving 

Ground  Smoke  tesT^  1977,  wISAa  Technical  Report  Th-2C)1,  Aberdeen  Proving 
Ground,  MD  (AD  A045874) 

^®F.  Pasquill,  1974,  Atmospheric  Diffusion,  second  edition,  Halsted  Press 
01 V.,  John  Wiley  and  Sons, Tnc.,  Hew  York 

^^C.  H.  B.  Priestley,  1956,  "A  Working  Theory  of  the  Bent-Over  Plume  of  Hot 
Gas,"  Quart  J  Roy  Meteorol  Soc,  82:165-176 
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The  vertical  (ambient)  temperature  profile  is  modeled  as 


T 


z 


1  -  exp[-a{z  -  Z-)/H  ] 

=  T  +  T„  - 1 

°  1  -  exp[-a(H^  -  2^)/HJ 


Z  <  H  /lO 
—  m 


T 


z 


=  Tio 


r ^hm  -  Tio  1 

I  Hm  -  ^ 


z  >  H  /lO 
m 


(30) 


Tio  =  VH^IO)  . 


where  T^.  is  the  temperature  measured  at  the  reference  height  (z^),  is  the 
temperature  at  the  mixing  height  (H^)i  and  a  is  chosen  so  as  to  fit  to  the 
measured  ambient  temperature  gradient  at  the  reference  height  and  to  the 
adiabatic  lapse  rate  (0.00966“C/m)  at  z  *  H^/10* 

The  vertical  velocity  (u)  at  any  later  time  is  found  by  first  applying  the 
conservation  of  momentum  along  the  vertical: 


d(u  g 

«  "  ■  '^m“  ’  (31) 

z 

where  AT  is  cloud  temperature  excess,  g  is  the  acceleration  due  to  gravity,  T^ 
as  before  is  the  (absolute)  ambient  air  temperature  at  the  centroid  height, 
and  is  the  momentum  mixing  coefficient  taken  to  be  0.10 

Equation  (31)  with  (29)  and  (30)  is  then  solved  for  co  using  reiterative  tech¬ 
niques  assuming  zero  initial  velocity  to  further  determine  the  rise  function 
by  way  of  the  following  kinematic  relations: 


u(t)  =  u(t  -  At)  +  (du)/dt)At; 
H(t)  =  H(t  -  At)  +  a)(t)At  . 


(32) 


H.  B.  Priestley,  1956,  "A  Working  Theory  of  the  Bent-Over  Plume  of  Hot 
Gas,"  Quart  J  Roy  Meteorol  Soc,  82:165-176 


The  process  is  then  repeated  by  incrementing  time  (and  hence  C)  over  the  "age" 
of  the  cloud.  Beyond  the  first  time  increment,  a  term  [(aT)  x  aC/C]  is 
sequentially  added  to  equation  (29)  to  account  for  the  vertically  changing 
temperature  of  the  ambient  entrained  air.  In  actual  practice  the  time  incre¬ 
ment  is  computed  so  as  to  limit  the  cloud  rise  increment  to  1  m  or  less  to 
assure  convergence  of  the  numerical  procedure. 


6.  EXAMPLE  (FROM  SMOKE  WEEK  II) 


6.1  Conversion  of  Model  Results 

As  mentioned  in  section  1,  the  model  was  coded  in  such  a  way  as  to  be  compat¬ 
ible  with  measurements  made  during  the  major  field  tests.  This  coding  allows 
nearly  direct  comparison  between  model  results  and  measured  data.  However,  a 
word  of  caution  is  required  to  interpret  the  results  appropriately. 

The  model  has  assumed  radiometric  units  throughout,  whereas  the  units  reported 
in  the  field  tests  are  mixed;  that  is,  sky  and  solar  data  are  in  radiometric 
units,  but  target,  background,  and  path  radiance  are  in  photopic  units. 
Because  the  underlying  spectrum  is  not  uniform  (that  is,  the  sky  is  blue, 
clouds  are  white,  and  the  sun  is  yellow-green),  some  error  and  confusion 
result  when  converting  between  the  two  systems.  Rather  than  try  to  correct 
for  the  nonuniform  spectrum  (a  procedure  which  could  only  increase  the  error), 
we  choose  here  to  assume  the  spectrum  nearly  uniform  and  convert  the  model 
results  to  photopic  units  by  using  the  standard  photopic  response  curve^'^ 
which  can  be  closely  approximated  by  the  following  Gaussian  function  (see 
figure  5): 


Rp  =  Rg  exp[-l/2(x  -  x^j)/o)2]  , 


(33) 


where 

Rq  =  673  Im/W  , 

x^  =  0.56um, 

0 

o  =  0.0426vim, 


22a.  Stimson,  1974,  Photometry  and  Radiometry  for  Engineers,  John  Wiley  and 
Sons,  Inc.,  New  York 


which  can  be  integrated  to  yield  the  following  conversion  factor: 


E(lm)  =  Rp/Tir  o  E(W) 


(34) 


or  in  terms  of  bandwidth  (full  width  at  half  maximum): 


E(lm)  =  Rq  ax  E(W)  . 


(35) 


Both  the  bandwidth  (ax)  and  position  of  maximum  response  (x^)  are  input  by  the 
user  and  are  O.lOOum  and  0.56um  for  straight  photopic  conversion.  For  input 
wavelengths  other  than  0.56um,  the  model  shifts  and  reduces  the  peak  response 
via  a  multiplicative  Gaussian  factor: 


f  =  exp 


(36) 


which  is  equal  to  unity  for  x  =  0.56  and  is  essentially  zero  for  the 
infrared.  In  all  cases  the  model  also  provides  output  in  radiometric  units. 

Also  the  Smoke  Week  sun  and  sky  radiances  are  reported  for  a  detector  field  of 
view  of  1®  requiring  division  by  (Tr/lSO)^  to  convert  to  a  unit  steradian. 
This  conversion  is  not  required  for  the  input  to  the  model  as  now  coded  since 
these  data  are  used  only  in  a  relative  sense.  For  sake  of  completeness,  some 
further  required  conversion  factors  are: 


1  footcandle  =  10.76  lumens  per  square  meter 


1  footlambert  =  candles  per  square  meter 


(37) 


1  candle  =  1  lumen  per  steradian 


6.2  Input  Data 

Trial  1  of  the  Smoke  Week  II  field  test,  held  at  Eglin  Air  Force  Base, 
Florida,  in  November  1978  consisted  of  the  detonation  of  15  155-mm  hexachlo- 
roethane  (HC)  Type  Ml  canisters  arranged  in  the  configuration  sketched  in 
figure  6.  The  source  characteristics  used  in  the  model  were  those  as  reported 
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in  EOSAEL  80  Technical  Documentation;^^  and  the  mass  extinction  coefficient, 
single  scattering  albedo,  and  phase  function  were  those  of  Shirkey,  Clayton, 
and  Quintis^®  for  HC  smoke.  For  modeling  purposes,  the  munitions  were 
separated  into  four  groups  as  indicated  by  the  sketched  outlines  in  figure  6 
with  each  group  treated  as  single-point  detonation  of  appropriate  total 
mass.  For  buoyant  smokes  this  latter  procedure  may  cause  some  concern; 
however,  for  HC  munitions  which  are  only  slightly  buoyant  this  procedure 
causes  only  insignificant  errors. 

Meteorological  conditions  during  the  test  were  typical  of  fair  weather  with  30 
percent  cloud  cover  and  11.3  km  visibility.  Model  inputs  either  taken  from 
the  original  test  report  or  derived  (estimated)  from  data  therein  are  listed 
in  table  2.  Table  3  lists  the  ambient  sky  radiation  measurements  made  during 
these  tests. 

The  sky  radiance  map  derived  from  the  data  of  table  3  for  the  model  sector 
midpoints  is  shown  in  figure  7. 

Figure  8  shows  the  modeled  and  measured  results  for  both  path  integrated 
concentration  (figure  8a)  and  path  luminance  (figure  8b).  The  results  for 
path  integrated  concentration,  although  overall  high,  are  typical  of  those 
reported  in  other  validation  studies. The  path  luminance  results  are  most 
interesting  in  that  the  brightening  effect  at  the  cloud  edges  is  quite  notice¬ 
able.  This  effect  is  often  observed  in  natural  clouds  and  is  referred  to  as  a 
"silver  lining."  The  occurrence  and  magnitude  of  the  bright  edges  depend 
strongly  upon  the  angular  distribution  of  ambient  radiation.  The  overall 
agreement  between  model  data  and  data  of  figure  8b  is  encouraging. 

Figure  9  is  a  more  detailed  examination  of  the  cloud  at  time  t  =  100  s.  Here 
both  the  direct  and  diffuse  components  of  radiation  are  plotted  as  a  function 
of  depth  of  penetration.  This  procedure  may  be  viewed  as  moving  the  target 
into  the  cloud  along  the  LOS  away  from  the  observer.  Until  a  significant 
port^"'!  of  the  cloud  is  penetrated,  the  diffuse  component  is  near  zero  and  the 


a.  Sutherland,  1981,  "Smoke  Obscuration  Model,"  chapter  3,  EOSAEL  80, 
Volume  1,  Technical  (locumentation,  editor  L.  D.  Duncan,  ASL-TR-0072,  uS  Army 
Atmospheric  Sciences  laboratory.  White  Sands  Missile  Range,  NM  (AD  B055130L) 

2®R.  C.  Shirkey,  D.  Clayton,  and  D.  M.  Quintis,  1981,  "Aerosol  Phase  Function 
Data  File  PFNDAT,"  chapter  16,  EOSAEL  80,  Volume  II,  Users  Manual,  editors 
R.  C.  Shirkey  and  S.  G.  O'Brien,  ASl-TR-0073,  US  Army  Atmospheric  Sciences 
Laboratory,  White  Sands  Missile  Range,  NM  (AD  B056119) 

2®D.  W.  Hoock,  R.  A.  Sutherland,  and  D.  Clayton,  1981,  Comparisons  Between  the 
EOSAEL  80  Model  SMOKE  and  the  Inventory  Munition  Test  r)^ase  lla.  Technical 
Report,  Us  Army  Atmospheric  Sciences  Laboratory,  White  Sands  Missile  Range,  NM 
(in  process) 


direct  component  Is  at  a  maximum.  As  the  target  moves  Into  the  cloud,  the 
diffuse  component  Increases  while  the  direct  component  decreases.  The  net 
effect  Is  a  reduction  In  both  the  direct  signal  and  the  contrast. 

A  discussion  of  how  these  (and  other)  outputs  of  the  model  can  be  used  In 
other  smoke  screening  and  perception  models  can  be  found  elsewhere.^° 


30d.  W.  Hoock  and  R.  A.  Sutherland,  1981,  "Path  to  Background  Luminance  Ratios 
for  the  EOSAEL  80  Munitions  Expenditure  Model  SCREEN,"  Proceedings  of  Smoke 
Symposium  V,  Harry  Diamond  Laboratories,  Adelphi,  MD 
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TABLE  1.  DIFFUSION  PARAMETERS  USED  IN  THE  TRANSPORT  AND  DIFFUSION  ROUTINE  FOR 
VARIOUS  VALUES  OF  ROUGHNESS  PARAMETER  (Z^)  AND  STABILITY  CATEGORY.* 


Default 


0  <  <  10 
—  0 


10  <  Zq  <  100 


Z„  >  100 

0 


STABILITY  CATEGORY 


A 

B 

C 

D 

E 

F 

Ax 

7.57 

7.57 

7.57 

7.57 

7.57 

7.57 

Bx 

0.93 

0.93 

0.93 

0.93 

0.93 

0.93 

Ay 

44.1 

44.1 

44.1 

24.8 

15.6 

15.6 

■“y 

1.50 

1.50 

1.50 

0.88 

0.66 

0.66 

Az 

16.4 

16.4 

16.4 

14.2 

12.7 

12.7 

Bz 

1.50 

1.50 

1.50 

0.88 

0.66 

0.66 

Ax 

2.77 

3.55 

5.38 

8.68 

12.6 

17.5 

Bx 

0.90 

0.90 

0.90 

0.90 

0.90 

0.90 

Ay 

2.77 

3.55 

5.38 

8.68 

12.6 

17.5 

0.90 

0.90 

0.90 

0.90 

0.90 

0.90 

»z 

7.32 

9.65 

12.0 

15.6 

26.3 

47.6 

Bz 

0.94 

0.89 

0.85 

0.81 

0.78 

0.72 

Ax 

2.77 

3.55 

5.38 

8.68 

12.6 

17.55 

Bx 

0.90 

0.90 

0.90 

0.90 

0.90 

0.90 

Ay 

2.77 

3.55 

5.38 

8.68 

12.6 

17.5 

By 

0.90 

0.90 

0.90 

0.90 

0.90 

0.90 

Az 

4.13 

5.78 

6.75 

8.59 

14.9 

24.6 

Bz 

0,90 

0,85 

0.80 

0.76 

0.73 

0.67 

Ax 

2.77 

3.55 

5.38 

8.68 

12.6 

17.5 

Bx 

0.90 

0.90 

0.90 

0.90 

0.90 

0.90 

N 

2.77 

3.55 

5.38 

8.68 

12.6 

17.5 

»y 

0.90 

0.90 

0.90 

0.90 

0.90 

0.90 

*z 

1.80 

2.23 

2.40 

3.17 

5.07 

7.57 

Bz 

0.83 

0.77 

0.72 

0.68 

0.65 

0.58 

*Zq  In  centimeters,  x,  y,  and  z  in  meters 


TABLE  2. 

AMBIENT  RADIATION  MEASUREMENTS 

FIELD  TEST.* 

(VISIBLE)  FROM  SMOKE  WEEK  II, 

TRIAL  1 

Sky  Radiance 

W  sr‘i 

m”2 

0  + 

3.7" 

70.7“ 

160.7“ 

250.7“ 

0“ 

6.2 

6.2 

6.2 

6.2 

10“ 

9.8 

12.1 

15.8 

13.1 

20“ 

7.5 

13.1 

20.7 

13.1 

30“ 

6.6 

5.1 

33.8 

13.1 

40“ 

7.5 

17.4 

54.5 

12,1 

50“ 

8,5 

19.7 

133.6 

13.1 

60“ 

10,8 

21.3 

124.7 

13.1 

O 

O 

13.1 

21.3 

124.7 

17.4 

o 

O 

CO 

17.4 

10.8 

103.4 

17.4 

90“ 

18.4 

8.9 

68.6 

18.1 

Solar: 

Zenith 

Azimuth 

Beam 

flux  (W/m^) 

51.5 

151.9 

208.7 

Surface  Irradiance  (W/m^):  571.5 


♦Note  that  reported  sky  radiance  data  must  be  divided  by  {n/180)2  to  convert 
unit  solid  angle  to  1  sr.  Also  a  factor  10"^  converts  uW/cm^  to  W/m^ 


TABLE  3.  METEOROLOGICAL  INPUTS  TO  THE  MODEL  FOR 
SMOKE  WEEK  II,  TRIAL  1  FIELD  TEST 


Windspeed  (8  m) 

Wind  direction  (8  m) 

Wind  power  law  exponent 

Ambient  air  temperature  (1.0  m) 

Temperature  gradient 

Mixing  height 

Mixing  height  temperature 

Stability  category  (Pasquill) 

Relative  humidity 

Dew-point  temperature 

Surface  irradiance  (short  wave) 

Surface  temperature 

Surface  reflectivity 

Surface  roughness 

Surface  particle  reflectance 


4.1  m/s 
116.3  deg 
0.11 
23.6*C 
-0.36“C/m 

400  m  (derived  in  model) 
18.3“C  (derived  in  model) 

C 

52% 

not  needed 
0.82  Langley /min 
24.0®C  (derived  in  model) 
0.25  (estimated  from  data) 
0.0  (default) 

1.0  (default) 
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Figure  1. 


Figure  2. 


Sketch  demonstrating  the  effects  of  extinction  and  path  radiance 
on  radiant  energy  received  by  an  observer. 

SOURCE  (i,  e.,  SUN)  — x 


Sketch  demonstrating  the  scattering  of  ambient  radiation  into 
the  LOS. 
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Figure  5.  Photopic  response  curve  and  Gaussian  functional  fit  for  converting 
radiance  (V/m^)  to  lualnance  (candles/a^). 
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ATTN:  R52/0r.  Douglas  Moods 
Fort  George  G.  Meade,  MD  20755 

Chief 

Intelligence  Materiel  Development 
&  Support  Office 

US  Army  Electronic  Warfare  Laboratory 
ATTN:  DELEW- I  (LTC  Kenneth  E.  Thomas) 
Fort  George  G.  Meade,  MD  20755 

The  Johns  Hopkins  University 
Applied  Physics  Laboratory 
ATTN:  Dr.  Michael  J.  Lun 
John  Hopkins  Road 
Laurel),  MO  20810 

Dr.  Stephen  T.  Hanley 
1720  Rhodesia  Avenue 
Oxon  Hill,  MD  20022 

Science  Applications  Inc. 

ATTN:  Mr.  6.  D.  Currie 
15  Research  Drive 
Ann  Arbor,  MI  48103 

Science  Applications  Inc. 

ATTN:  Dr.  Robert  E.  Turner 
15  Research  Drive 
Ann  Arbor,  MI  48103 

Commander 

US  Army  Tank -Automotive  Research 
S  Development  Command 
ATTN:  DRDTA-ZSC  (Mr.  Harry  Young) 
Warren,  MI  48090 

Commander 

US  Army  Tank  Automotive  Research 
A  Development  Command 
ATTN:  DRDTA-ZSC  (Mr.  Wallace  Mick,  Jr.) 
Warren,  MI  48090 


Dr.  A.  D.  Belmont 
Research  Division 
Control  Daca  Corporation 
P.O.  Box  1249 
Minneapolis,  MN  55440 

Director 

US  Army  Engr  Waterways  Experiment  Station 
ATTN:  MESEN  (Mr.  James  Mason) 

P.O.  Box  631 
Vicksburg,  MS  39180 

Dr.  Jerry  Davis 
Department  of  Marine,  Earth 
and  Atmospheric  Sciences 
North  Carolina  State  University 
Raleigh,  NC  27650 

Commander 

US  Arny  Research  Office 
ATTN:  DRXRO-GS  (Dr.  Leo  Alpert) 

P.O.  Box  12211 

Research  Triangle  Park,  NC  27709 
Commander 

US  Arnv  Research  Office 
ATTN:  ORXRO-PP  (Brenda  Mann) 

P.O.  Box  12211 

Research  Triangle  Park,  NC  27709 
Commander 

US  Arirty  Cold  Regions  Research 
A  Engineering  Laboratory 
ATTN:  CRREL-RD  (Dr.  K.  F.  Sterrett) 
Hanover,  NH  03755 

Commander/Director 
US  Amy  Cold  Regions  Research 
*  Engineering  Laboratory 
ATTN:  CRREL-RG  (Mr.  George  Aitken) 
Hanover,  NH  03755 

Commander 

US  Amy  Cold  Regions  Research 
S  Engineering  Laboratory 
AHN:  CRREL-RG  (Mr.  Roger  H.  Berger) 
Hanover,  NH  03755 

Commander 

US  Amy  Armament  Research 
i  Development  Command 
ATTN:  DRDAR-AC  (Mr.  James  Greenfield) 
Dover,  NJ  07801 
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Commander 

US  Am^  Armament  Research 
&  Development  Command 
ATTN:  OROAR-TSS  (Bldg  #59) 

Dover,  NJ  07801 

Commander 

US  Army  Armament  Research 
&  Development  Command 
AHN:  DRCPM-CAWS-EI  (Mr.  Peterls  Jansons) 
Dover,  NJ  07801 

Commander 

US  Army  Armament  Research 
A  Development  Command 
ATTN:  DRCPM-CAWS-EI  (Mr.  G.  H.  Waldron) 
Dover,  NJ  07801 

Deputy  Joint  Project  Manager 
for  Navy/USMC  SAL  GP 

ATTN:  DRCPM-CAWS-NV  (CPT  Joseph  Mi  cell) 
Dover,  NJ  07801 

Cotnmander/Di  rector 

US  Army  Combat  Surveillance  S  Target 
Acquisition  Laboratory 
ATTN:  DELCS-I  (Mr.  David  Longinotti) 

Fort  Monmouth,  NJ  07703 

Commander/Di  rector 

US  Army  Combat  Surveillance  4  Target 
Acquisition  Laboratory 
ATTN:  DELCS-PE  (Mr.  Ben  A.  Di  Campli) 

Fort  Monmouth,  NJ  07703 

Commander/Di rector 

US  Army  Combat  Surveillance  4  Target 
Acquisition  Laboratory 
ATTN:  DELCS-R-S  (Mr.  Donald  L.  Foiani) 
Fort  Monmouth,  NJ  07703 

Di rector 

US  Army  Electronics  Technology  4 
Devices  Laboratory 
AHN:  DELET-DD  (S.  Danko) 

Fort  Monmouth,  NJ  07703 

Project  Manager 
FIREFINOER/REMBASS 

ATTN:  DRCPM-FFR-TM  (Mr.  John  M.  BlaVo) 
Fort  Monmouth,  NJ  07703 


Commander 

US  Amy  Electronics  Research 
4  Development  Command 
AHN:  DRDEL-SA  (Dr.  Walter  S.  McAfee) 
Fort  Monmouth,  NJ  07703 

OLA,  2WS  (MAC) 

Holloman  AFB,  NM  88330 
Commander 

Air  Force  Weapons  Laboratory 
AHN:  AFWL/WE  (MW  John  R.  El  rick) 
Kirtland,  AFB,  NM  87117 

Di rector 

USA  TRADOC  Systems  Analysis  Activity 
ATTN:  ATAA-SL 

White  Sands  Missile  Range,  NM  88002 
Di rector 

USA  TRADOC  Systems  Analysis  Activity 
ATTN:  ATAA-SL  (Dolores  Anguiano) 

White  Sands  Missile  Range,  88002 

Di rector 

USA  TRADOC  Systems  Analysis  Activity 
ATTN:  ATAA-TDB  (Mr.  Louie  Dominguez) 
White  Sands  Missile  Range,  NM  88002 

Director 

USA  TRADOC  Systems  Analysis  Activity 
ATTN:  ATAA-TDB  (Mr.  William  J.  Leach) 
White  Sands  Missile  Range,  NM  88002 

Director 

USA  TRADOC  Systems  Analysis  Activity 
AHN:  ATAA-TGP  (Mr.  Roger  F.  Willis) 
White  Sands  Missile  Range,  NM  88002 

Director 

Office  of  Missile  Electronic  Warfare 
AHN:  DELEW-M-STO  (Dr.  Steven  Kovel) 
White  Sands  Missile  Range,  NM  88002 

Office  of  the  Test  Director 
Joint  Services  EO  GW  CM  Test  Program 
AHN:  DRXDE-TD  (Mr.  Weldon  Findley) 
White  Sands  Missile  Range,  NM  88002 

Commander 

US  Amy  White  Sands  Missile  Range 
AHN:  STEWS-PT-AL  (Laurel  B.  Saunders) 
White  Sands  Missile  Range,  NM  88002 
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Conmander 

US  Am^  RSD  Coordinator 
US  Embassy  -  Bonn 
Box  165 

APO  New  York  09080 

Grumman  Aerospace  Corporation 
Research  Department  -  MS  A08-35 
ATTN:  John  E.  A.  Selby 
Bethpage,  NY  11714 

Rome  Air  Development  Center 
ATTN:  Documents  Library 
TSLD  (Bette  Smith) 

Griffiss  AFB,  NY  13441 

Or.  Roberto  Vaglio-Laurin 
Faculty  of  Arts  and  Science 
Dept,  of  Applied  Science 
26-36  Stuyvesant  Street 
New  York.  NY  10003 

Air  Force  Wright  Aeronautical  Laboratories/ 
Avionics  Laboratory 

ATTN'  AFWAL/AARI-3  (Mr.  Harold  Geltmacher) 
Wright-Patterson  AFB,  OH  45433 

Air  Force  Wright  Aeronautical  Laboratories/ 
Avionics  Laboratory 

ATTN:  AFWAL/AARI-3  (CPT  William  C.  Smith) 
Wright-Patterson  AFB,  OH  45433 

Commandant 

US  Army  Field  Artillery  School 
ATTN:  ATSF-CF-R  (CPT  James  M.  Watson) 

Fort  Sill,  OK  73503 

Commandant 

US  Army  Field  Artillery  School 
ATTN:  ATSF-CD-MS 
Fort  Sill,  OK  73503 

Commandant 

US  Army  Field  Artillery  School 
ATTN:  ATSF-CF-R 
Fort  Sill,  OK  73503 

Commandant 

US  Army  Field  Artillery  School 
ATTN:  NOAA  Liaison  Officer 

(CDR  Jeffrey  G.  Carlen) 

Fort  Sill,  OK  73503 


Commandant 

US  Army  Field  Artillery  School 
Morris  Swett  Library 
ATTN:  Reference  Librarian 
Fort  Sill,  OK  73503 

Commander 

Naval  Air  Development  Center 
ATTN:  Code  301  (Mr.  George  F.  Eck) 
Warminster,  PA  18974 

The  Univei'-sity  of  Texas  at  El  Paso 
Electrical  Engineering  Department 
ATTN:  Dr.  Joseph  H.  Pierluissi 
El  Paso,  7X  79968 

Commandant 

US  Army  Air  Defense  School 

AHN:  ATSA-CD-SC-A  (CPT  Charles  T.  Thorn) 

Fort  Bliss.  TX  79916 

Commander 

HQ,  TRADOC  Combined  Arms  Test  Activity 
ATTN:  ATCAT-OP-Q  (CPT  Henry  C.  Cobb,  Jr.) 
Fort  Hood,  TX  76544 

Commander 

HQ,  TRADOC  Combined  Arms  Test  Activity 
AHN;  ATCAT-SCI  (Dr.  Darrell  W.  Collier) 
Fort  Hood,  TX  76544 

Commander 

US  Army  Dugway  Proving  Ground 
ATTN:  STEDP-MT-DA-L 
Dugway ,  UT  84022 

Commander 

US  Army  Dugway  Proving  Ground 

ATTN:  STEDP-MT-DA-M  (Mr.  Paul  E.  Carlson) 

Dugway,  UT  84022 

Commander 

US  Arny  Dugw4y  Proving  Ground 

ATTN:  STEDP-MT-DA-T  (Mr.  John  Trethewey) 

Dugway ,  UT  84022 

Commander 

US  Amy  Dugwuy  Proving  Ground 

AHN:  STEDP-MT-OA-T  (Mr.  William  Peterson) 

Dugway,  UT  84022 
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Defense  Documentation  Center 
ATTN:  ODC-TCA 
Cameron  Station  Bldg  5 
Alexandria,  VA  22314 
12 

Ballistic  Missile  Defense  Program  Office 
ATTN:  DACS-BMT  (Colonel  Harry  F.  Ennis) 
5001  Eisenhower  Avenue 
Alexandria,  VA  22333 

Defense  Technical  Information  Center 
ATTN:  DDA-2  (Mr.  James  E.  Shafer) 
Cameron  Station,  Bldg  5 
Alexandria,  VA  22314 

Commander 

US  Army  Materiel  Development 
&  Readiness  Command 

ATTN:  DRCBSI-EE  (Mr.  Albert  Giambalvo) 
5001  Eisenhower  Avenue 
Alexandria,  VA  22333 

Commander 

US  Army  Materiel  Development 
4  Readiness  Command 
ATTN:  DRCLOC  (Mr.  James  Bender) 

5001  Eisenhower  Avenue 
Alexandria,  VA  22333 

Defense  Advanced  Rsch  Projects  Agency 
ATTN:  Steve  Zakanyez 
1400  Wilson  Blvd 
Arlington,  VA  22209 

Defense  Advanced  Rsch  Projects  Agency 
ATTN:  Dr.  James  Tegnelia 
1400  Wilson  Blvd 
Arlington,  VA  22209 

Institute  for  Defense  Analyses 
ATTN:  Mr.  Lucien  M.  Biberman 
400  Army-Navy  Drive 
Arlington,  VA  22202 

Institute  for  Defense  Analyses 
ATTN:  Dr.  Ernest  Bauer 
400  Army-Navy  Drive 
Arlington,  VA  22202 

Institute  for  Defense  Analyses 
ATTN:  Dr.  Hans  G.  Wolf hard 
400  Army-Navy  Drive 
Arlington,  VA  22202 


System  Planning  Corporation 
ATTN:  Mr.  Daniel  '^riedman 
1500  Wilson  Boulevard 
Arlington,  VA  22209 

System  Planning  Corporation 
ATTN;  COL  Hank  Shelton 
1500  Wilson  Boulevard 
Arlington,  VA  22209 

US  Army  Intelligence  *  Security  Command 
ATTN:  Edwin  Speakman,  Scientific  Advisor 
Arlington  Hall  Station 
Arlington,  VA  22212 

Commander 

US  Army  Operational  Test 
&  Evaluation  Agency 
ATTN;  CSTE-ED  (Mr.  Floyd  I.  Hill) 

5600  Columbia  Pike 
Falls  Church,  VA  22041 

Commander  and  Director 
US  Army  Engineer  Topographic  Laboratories 
ATTN:  ETL-GS-A  (Mr.  Thomas  Neidringhaus) 
Fort  Bel  voir,  VA  22060 

Di rector 

US  Army  Night  Vision  & 

Electro-Optics  Laboratory 
ATTN:  DELNV-L  (Dr,  Rudolf  G.  Buser) 

Fort  Bel  voir,  VA  22060 

Director 

US  Army  Night  Vision  & 

Electro-Optics  Laboratory 
ATTN:  DELNV-L  (Dr.  Robert  S.  Rodhe) 

Fort  Bel  voir,  VA  22060 

Director 

US  Army  Night  Vision  4 
Electro-Optics  Laboratory 
ATTN:  DELNV-VI  (Mr.  Joseph  R.  Moulton) 
Fort  Bel  voir,  VA  22060 

Di rector 

US  Army  Night  Vision  4 
Electro-Optics  Laboratory 
ATTN:  DELNV-VI  (Luanne  P.  Obert) 

Fort  Bel  voir,  VA  22060 


Director 

US  Army  Night  Vision 
A  Electro-Optics  Laboratory 
ATTN:  DELNV-VI  (Mr.  Thomas  W.  Cassidy) 
Fort  Bel  voir,  VA  22060 

Di rector 

US  Army  Night  Vision  & 

Electro-Optics  Laboratory 
ATTN:  DELNV-VI  (Mr.  Richard  J.  Bergemann) 
Fort  Bel  voir,  VA  22060 

Di rector 

US  Army  Night  Vision  & 

Electro-Optics  Laboratory 
ATTN;  DELNV-VI  (Dr.  James  A.  Ratches) 

Fort  Bel  voir,  VA  22060 

Commander 

US  Army  Training  &  Doctrine  Command 

ATTN:  ATCD-AN 

Fort  Monroe,  VA  23651 

Commander 

US  Army  Training  A  Doctrine  Command 
ATTN;  ATCD-AN -M 
Fort  Monroe,  VA  23651 

Commander 

US  Army  Training  &  Doctrine  Command 
ATTN:  ATCD-F-A  (Mr.  Chris  O'Connor,  Jr.) 
Fort  Monroe,  VA  23651 

Commander 

US  Army  Training  i  Doctrine  Command 
ATTN;  ATCO-IE-R  (Mr.  David  M.  Ingram) 
Fort  Monroe,  VA  23651 

Commander 

US  Army  Training  4  Doctrine  Command 
ATTN;  ATCD-M-I/ATCD-M-A 
Fort  Monroe,  VA  23651 

Commander 

US  Army  Training  &  Doctrine  Command 
ATTN;  ATDOC-TA  (Dr.  Marvin  P.  Pastel) 
Fort  Monroe,  VA  23651 

Department  of  the  Air  Force 
OL-I,  AWS 

Fort  Monroe,  VA  23651 


Department  of  the  Air  Force 
HQS  5  Weather  Wing  (MAC) 

ATTN:  5  WW/DN 

Langley  Air  Force  Base,  VA  23655 
Commander 

US  Army  JNSCOM/Quest  Research  Corporation 
ATTN:  Mr.  Donald  Wilmot 
6845  Elm  Street,  Suite  407 
McLean,  VA  22101 

General  Research  Corporation 
ATTN;  Dr.  Ralph  Zirkind 
7655  Old  Springhouse  Road 
McLean,  VA  22102 

Science  Applications,  Inc. 

8400  Westpark  Drive 
ATTN;  Dr.  John  E.  Cockayne 
McLean,  VA  22102 

US  Army  Nuclear  4  Chemical  Agency 
ATTN;  MONA-WE  (Dr.  John  A.  Berberet) 

7500  Backlick  Road,  Bldg  2073 
Springfield,  VA  22150 

Director 

US  Army  Signals  Warfare  Laboratory 
ATTN:  DELSW-EA  (Mr.  Douglas  Harkleroad) 
Vint  Hill  Farms  Station 
Warrenton,  VA  22186 

Di rector 

US  Army  Signals  Warfare  Laboratory 
ATTN:  DELSW-OS  (Dr.  Royal  H.  Burkhardt) 
Vint  Hill  Farms  Station 
Warrenton,  VA  22136 

Commander 

US  Army  Cold  Regions  Test  Center 
ATTN:  STECR-TD  (Mr.  Jerold  Barger) 

APO  Seattle,  WA  98733 

HQDA  (SAUS-OR/Hunter  M.  Woodall,  Jr./ 

Dr.  Herbert  K.  Fallin) 

Rm  2E  614,  Pentagon 
Washington,  DC  20301 

COL  Elbert  W.  Friday,  Jr. 

OUSORE 

Rm  3D  129,  Pentagon 
Washingtotr;  DC  20301 
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Defense  Communications  Agency 
Technical  Library  Center 
Code  222 

Washington,  DC  20305 
Di rector 

Defense  Nuclear  Agency 

ATTN:  Technical  Library  (Mrs.  Betty  Fox) 

Washington,  DC  20305 

Di rector 

Defense  Nuclear  Agency 
ATTN:  RAAE  (Or.  Carl  Fitz) 

Washington,  DC  20305 

Di rector 

Defense  Nuclear  Agency 

ATTN:  SPAS  (Mr.  Donald  J.  Kohler) 

Washington,  DC  20305 

Defense  Intelligence  Agency 
ATTN:  DT/AC  (LTC  Robert  Poplawski) 
Washington,  DC  20301 

HQDA  (DAf^-ARZ-O/Or.  Verderame) 
Washington,  DC  20310 

HQDA  (DAMI-ISP/Mr.  Beck) 

Washington,  DC  20310 

Department  of  the  Army 
Deputy  Chief  of  Staff  for 
Operations  and  Plans 
ATTN:  DAMO-RQ 
Washington,  D(i  20310 

Department  of  the  Army 
Director  of  Telecommunications  and 
Command  and  Control 
ATTN:  DAMO-TCZ 
Washington,  DC  20310 

Department  of  the  Army 

Assistant  Chief  of  Staff  for  Intelligence 

ATTN:  DAMI-TS 

Washington,  DC  20310 

HQDA  (DAEN-ROM/Dr.  de  Percin) 

Casimir  Pulaski  Building 
20  Massachusetts  Avenue 
Room  6203 

Washington,  DC  20314 
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National  Science  Foundation 
Division  of  Atmospheric  Sciences 
ATTN:  Or.  Eugene  W.  Bierly 
1800  G.  Street.  N.W. 

Washington,  DC  20550 

Director 

Naval  Research  Laboratory 

ATTN:  Code  4320  (Or.  Lothar  H.  Ruhnke) 

Washington,  DC  20375 

Commanding  Officer 

Naval  Research  Laboratory 

ATTN;  Code  6009  (Or.  John  MacCallum,  Jr.) 

Washington,  DC  20375 

Commanding  Officer 

Naval  Research  Laboratory 

ATTN;  Code  6530  (Mr.  Raymond  A.  Patten) 

Washington,  DC  20375 

Commanding  Officer 

Naval  Research  Laboratory 

ATTN;  Code  6533  (Dr.  James  A.  Dowling) 

Washington,  DC  20375 
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